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Abstract 

An epitaxial transition from a bicontinious double gyroid to a hexagonally packed cylinder 
structure induced by an external flow is simulated using real-space dynamical self-consistent 
field technique. In order to simulate the structural change correctly, we introduce a system 
size optimization technique by which emergence of artificial intermediate structures are sup- 
pressed. When a shear flow in [111] direction of the gyroid unit cell is imposed, a nucleation 
and growth of the cylinder domains is observed. We confirm that the generated cylindrical 
domains grow epitaxially to the original gyroid domains as gyroid <i{220} "~ > cylinder <i{io}. 
In a steady state under the shear flow, the gyroid shows different reconnection processes 
depending on the direction of the velocity gradient of the shear flow. A kinetic pathway pre- 
viously predicted using the self-consistent field theory where three fold junctions transform 
into five fold junctions as an intermediate state is not observed. 

KEYWORDS: diblock copolymer, bicontinuous double gyroid, self-consisten field theory, 
dynamical mean field theory, epitaxial transition 

1 Introduction 

In the past, the self-organized microdomain structures of diblock copolymers have been the 
target of extensive studies El El El 111 Especially, the order-order transitions (OOTs) between 
the microdomain structures are one of the central issues of the current experimental and 
theoretical studies. Among various microdomain structures, the bicontinuous double gyroid 
(G) structure has attracted a great interest because of its complex structure (space group 
Ia3d) 5 . Although this G phase exists only in a narrow region of the phase diagram of 
a diblock copolymer, i.e. the region between the lamellar (L) phase and the hexagonally 
packed cylinder (C) phase El, its complex domains are expected to have a wide applicability 
to various techniques, for example, microporous systems, nano-reactors, and so on El El El 

Experimentally, the OOT from the G phase to the C phase is believed to be an epitaxial 
transition where the created cylindrical domains are commensurate to the original gyroid 
domains. However, the microscopic detailed process of this epitaxial transition has not been 
understood yet. 
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Epitaxial relationships between the G and C structures were observed in several exper- 
iments. Upon a temperature change, Rangon and Charvolin have observed that the {10} 
plane of the C is in commensurate to the {211} plane of the G domains in a surfactant 
system^. In the present paper, we will use simple notation as C {10} — > G {211} for such 
an epitaxial relationships between planes. 

An external shear flow also accelerates the OOTs. Under a shear flow and a temperature 
change, Schulz et al. have found different epitaxial relationships C {10} — ► G {220} and C 
{11} — >• G {211} in a block copolymer mixture by small-angle neutron scattering (SANS) 
experiments ™ Under similar experimental conditions using small-angle X-ray scattering 
(SAXS) diffraction techniques, Forster et al. have observed the epitaxial relationship C 
{10} — > G {211} similar to the Rangon and Charvolin's observation^. The same epitaxial 
relationship was also observed by Vigild et al. in a block copolymer system using SANS 
A cyclic transition C — > G — ► C in a block copolymer solution has been studied by Wang 
and Lodge, who supported the epitaxial relationship G {211} — ► C {10}^. 

From the point of view of kinetics, a long-lived coexistence between the C phase and the 
G phase has been found in the C — > G of a block copolymer system under a shear flow and 
a temperature change Furthermore, a grain boundary between the C phase and the G 
phase has been observed by a polarized optical microscopy in a quenched polymer solution 
^3. These observations suggest an existence of a stable boundary between the C phase and 
the G phase. 

On the theoretical side, mean field theories have been used to investigate microdomain 
structures of diblock copolymers. Using the self-consistent field (SCF) technique, Helfand 
and Wasserman have evaluated the free energy and predicted the equilibrium domain sizes of 
the classical phases in the strong segregation regime such as the body centered cubic crystal 
of spherical domains (BCC), the C phase, and the lamellar (L) phase EHS3ED On the other 
hand, the phase diagram of diblock copolymer in the weak segregation regime was predicted 
by Leibler using the random phase approximation (RPA) ESI Leibler's phase diagram is com- 
posed of classical phases and the disordered (D) phase depending on the values of the block 
ratio and the '■> i- e - the product of the Flory-Huggins interaction parameter \ and the 
total degree of polymerization of diblock copolymer N. The entire phase diagram including 
both the weak segregation regime and the strong segregation regime has been constructed by 
Matsen and Shick using the SCF technique in the reciprocal lattice space. Besides the clas- 
sical phases, they predicted the complex G phase in the weak and intermediate segregation 
regime El ED This theoretical phase diagram was confirmed experimentally E2 

Despite the success of the mean field theories on the equilibrium phase behavior, the 
investigation on the dynamic properties has not been fully developed yet. There have been 
a few trials on the dynamics of OOTs and order-disorder transitions (ODTs) of the mi- 
crodomain structures of diblock copolymers using the mean field approximation. A time 
dependent Ginzburg-Landau (TDGL) model was used to investigate the instability in the 
OOTs and ODTs such as OOTs L -> C and L -> S, and C -> S. ESI El 123 In these studies, 
the authors retained the most unstable modes in the Fourier amplitudes of the density fluc- 
tuations emerging in the vicinity of the critical point. A TDGL model described in terms 
of Fourier modes with two sets of wave vectors with different magnitudes has been used to 
study the transitions D — ► G, G — > C, and so on ESI EH ESI Although the TDGL theory is 
efficient in investigating large-scale systems, it is in principle applicable only to the weak 
segregation regime. 

On the other hand, the SCF theory can be used to study the phase transitions in weak, 



intermediate and strong segregation regimes. The quantitative accuracy of the SCF theory is 
another advantage compared to the TDGL theory. This is because the SCF theory takes the 
conformational entropy of the polymer chains into account precisely E3 ESS ESS ED Using the 
SCF theory, Laradji et al. have investigated the epitaxial transitions such as L *-* C, C <-> S, 
and G — > C taking the anisotropic fluctuations into account E3E21 Matsen has also studied 
the transitions C <-> S and C G using the SCF theory and has proposed a nucleation and 
growth model of the epitaxial transitions^"^. 

All of these theoretical studies mentioned above rely on the reciprocal space representa- 
tions EHIE21ESI anc [ m ost of experimental studies EflEHED h ave supported the existence of the 
epitaxial OOT G {211} <-» C {10} except for Schulz et al., who have supported the epitaxial 
OOT G {220} «-> C {10}. Experimentally, the epitaxial OOT G {211} C {10} and G 
{220} <-> C {10} are recognized as the same epitaxial relationships ^because the diffraction 
peaks from both G {211} and G {220} match well with the diffraction peaks from the C 
{10}. In this argument, however, the kinetic pathway in real space was not considered. In 
Figure 1, we show a projection of the G structure onto the [111] direction in the real space. 
The epitaxial relations G {211} <-> C {10} and G {220} «-> C {10} are shown in Figure 
1(a) and 1(b), respectively, where the spacing of the G planes and the epitaxial cylindrical 
domains are shown. As the directions of these two planes G {112} and G {220} are perpen- 
dicular with each other and the spacings between adjacent planes are also different, the two 
growth mechanisms shown in Figures 1(a) and (b) should be regarded as different ones. 

Furthermore, the OOTs of block copolymer melts are first order phase transitions and 
the nucleation and growth process of domains is expected. Since these process is spatially 
inhomogeneous, such a transition is not compatible with the treatment in the reciprocal lat- 
tice space (Fourier space), where spatially periodic lattice structures are assumed. Therefore, 
dynamical simulations in real space, such as the dynamical SCF simulation, is necessary to 
correctly investigate this transition E3EZIEHIE2I 

In the present paper, we study the epitaxial OOT G — > C using the dynamical SCF 
theory under shear flows. In order to treat the first order transition, we introduce a system 
size optimization (SSO) method, in which the side lengths of the simulation box are auto- 
matically adjusted so that the size and the shape of the simulation box are fit for the lattice 
spacing and the lattice axes of the ordered structures. Recently, Barrat et al. proposed a 
similar technique to study equilibrium domain morphology of block copolymer systems 
In the course of the transition, we observe a complex transient state that is composed of 
cylindrical domains parallel to the G {220} plane. We also confirm that our SSO method 
can reproduce spatially inhomogeneous nucleation and growth processes. Actually we ob- 
serve the coexistence between the G phase and the C phase, which is consistent with the 
experimentally observed first order transition behavior^. Furthermore we found that the G 
structure shows different deformation behaviours depending on the direction of the velocity 
gradient of the shear flow. 

We cannot confirm the scenario of the transition proposed by Matsen^, where the three 
fold junctions transform into five fold junctions. 

Finally, we clarify the kinetic pathway from the G phase to the C phase under a shear 
flow. To the best of our knowledge, this kinetic pathway in the real space has not been 
reported in the literature. 



2 Theory 



2.1 Dynamical self-consistent field theory 

Here, we briefly summarize the SCF theory for an A-B diblock copolymer EHHHHHIIII Le^ 
us consider a melt of A-B diblock copolymer. Due to the screening effect in the melts, we 
can assume Gaussian statistics for the chain conformation. Within this Gaussian statistics, 
the K-tjpe (K =A or B) segment is characterized by the effective bond length bjc, and the 
.fT-type block is characterized by the degree of polymerization Nk- Then the total degree 
of polymerization N is defined as iV = Na + Nb- We introduce an index s to specify each 
segment, where s = corresponds to the free end of the A-block and s = N corresponds to 
the other free end of the I?-block. Therefore, < s < Na and iV^ < s < N correspond to the 
A- block and the 5-block, respectively. In order to evaluate the conformational entropy, we 
need the statistical weight of any subchains. Let us use the notation Q(s' , r'; s, r) to denote 
the statistical weight of a subchain between s-th and s'-th segments (0 < s' < s < N) that 
are fixed at the positions r and r'. This statistical weight can be obtained by solving the 
following Edwards equation within the mean-field approximation 

^-Q(s',r';s,r)= \ b -^V 2 - (3V(s,r)]Q(s' ,r'; s,r), (1) 



ds 
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where j3 = 1/(/cbT), b(s) = bx if the s-th segment is the K-type segment, and V(s,r) is an 
external potential acting on the s-th segment at r imposed by the surrounding segments. 
Here, we assume that the external potential V(s, r) is the same if the segment species (A or 
B) is the same. Thus, 

, \ _ J Va(t) if s indicates an A-segment . . 

1 Vg(r) if s indicates a 5-segment. 

Equation (0) should be supplemented by the initial condition Q(0, r'; 0, r) = <5(r' — r). As the 
two ends of the block copolymer are not equivalent, we should introduce another statistical 
weight Q(s', r'; s, r), which is calculated in the opposite direction along the chain starting 
from the free end s = N. 

To reduce the computational cost, we define an integrated statistical weights q(s, r) and 
q(s, r)as follows: 

q{s,r) = J dr'Q(0,r'; S ,r) 

q(s,r)= fdr'Q(0,r';s,r). (3) 



It is easy to confirm that q(s, r) and q(s, r) also satisfy eq. (P). 

By using eqs. ©, the density of the K-type segments at position r is given by 

<p K (r) = C f ds q(s,r)q{N - s,r), (4) 

J s£K -block 



where C is the normalization constant: 



The parameter M is the total number of chains in the system and Z is the single chain 
partition function which is independent of K, i.e. Z — j drq(s, r)q(N — s, r) = / drq(N, r) = 
Jdrq(N,T). 

The external potential Vk{t) can be decomposed into two terms as follows 

V K (r) = £KK><t>K>{?) ~ /^(r). (6) 

K' 

The first term is the interaction energy between segments and the €kk' is the nearest- 
neighbor pair interaction energy between a i^-type segment and a K'-type segment, which 
is related to the Flory-Huggins interaction parameter via xab = z /3[^ab — (1/2) (e^A + £bb)] 
where z is the number of nearest neighbor sites. The /iR-(r) is the chemical potential of 
the if- type segment, which is the Lagrange multiplier that fixes the density of the K-type 
segments at the position r to the specified density value. The Vk(t) must be determined in 
a self-consistent manner so that this constraint is satisfied. Such a self-consistent condition 
is achieved by an iterative refinement of the Vjf(r). 

To improve the stability of the numerical scheme, we used the following finite difference 
scheme for the Edwards equation, eq. (JIJ 

r 0V(s,r)As,/^ b(s) 2 „ 9A \ r 8V(s.t)As, , s 
q(s + As, r) = exp[- ^ 1 1 ; ] (l + -±J_V 2 As) exp[- P V ^ ) }q(s, r). (7) 

The Helmholtz free energy of the system can be given as follows 

r=-k B TM\nZ + \Y.H I ' dT6 KK ,<p K {r)<p K ,{v) - £ f dvV K {v)<p K {r). (8) 

K K 1 K J 

To introduce dynamics into the model, we assume Fick's law of linear diffusion for the 
segment densities and an effect of the flow advection as follows 

J^<Mr, t) = L K VVW - V{v(r, t)<j> K {r, *)}, (9) 

where Lk is the mobility of K-type segment and v(r, t) is the local flow velocity such as the 
velocity of the externally imposed shear flow. 



2.2 System size optimization method 

Periodic microdomain structures of diblock copolymers have the crystal symmetry. To ob- 
tain equilibrium states of these periodic structures using the mean field theory, the free 
energy density of the system must be minimized with respect to the lattice structures of 
the ordered microdomains. Same is true for two phase coexisting states where the system 
size should be optimized with respect to the coexisting two periodic structures. For these 
purpose we introduce the system size optimization (SSO) method that minimizes the free 
energy density of the system by optimizing the side lengths of the simulation box on which 
periodic boundary conditions are imposed. This is a similar method as the constant pressure 
molecular dynamics simulation proposed by AnderserP^. In the static SCF calculations, this 
optimization can be performed by requiring the following local equilibrium condition for each 
side length of the simulation box: 

I- - < 10 > 



where Ci(i = x,y, z) is the side length of the simulation box. The left-hand side of eq. (fTU|) 
can be evaluated numerically using the following central difference approximation 



dT _ FjCj + ACj) - FjCj - AC,) 

where ACi is a small variation of Li. We used the parabolic optimization method^ to solve 
eq. (HUJ). 

On the other hand, when the dynamical SCF calculation is performed, we should regard 
Ci as a dynamical variable whose dynamics is described by the following ficticious equation 
of motion 

dCi dT 

~dt = ~ Cl ac? (12) 

where Q is a positive coefficient whose value is chosen properly so that the local equilibrium 
condition eq-lfTUj) for is guaranteed at every time step. 

We checked the validity of our dynamical SSO method by using an A-B diblock copoly- 
mer melt whose stable equilibrium phase is the C phase. We performed two dimensional 
simulations where we assumed that ( x = ( y = ( for simplicity, and we changed £ from 0.0 
to 0.5. The parameters characterizing the A-B diblock copolymer are as follows: the total 
length of the copolymer N = 20, the block ratio of the A block / = Na/N = 0.35, and the 
effective bond lengths of each segment type are unity. The interaction parameter is set to 
be \N = 15, which corresponds to the C phase in its equilibrium state ™ The initial state 
is set to the D phase to which we added small random noise with the standard deviation 
0.0006. The initial shape of the simulation box is a square with side length 32.0. As the 
square shape of the simulation box is not compatible with the perfect C phase, the SSO 
method adjusts the side lengths of the simulation box automatically. 

In Figure 2, we show a comparison of the domain morphologies in the late stage (t = 5000) 
between the two cases (a) with ( = 0.001 and (b) with ( = 0.05, respectively. In case (a), the 
C structure is distorted because the rate of the change in the side lengths of the simulation 
box is too slow to catch up with the change in the domain periodicity. On the other hand 
in case (b), a perfect C phase is realized. When £ = 0.5, we observed that the dynamical 
scheme eq. (|T2"|) becomes unstable. 

Other dynamical variables that depend on the value of ( are shown in Figures 3 and 4. 
Figure 3 shows the time evolution of the free energy. The dotted line is the reference state 
with C = (i.e. the case without SSO) which reaches the distorted morphology shown in 
Figure 2(a). When the value of ( is small (( = and 0.001), it takes longer time for the free 
energy to relax and finally the system is trapped in a local minimum of the free energy. For 
the intermediate values of ( (( = 0.05, 0.1, and 0.2), the system reaches the perfect C phase 
as shown in Figure 2(b). When the value of ( is large (( = 0.5), the free energy initially 
drops rapidly and then the system is trapped by a local minimum of the free energy. These 
results mean that choosing an appropriate value of ( accelerates the system to relax to the 
equilibrium domain morphology without distortions and defects. 

Figure 4 shows the time evolution of the side lengths of the simulation box. The solid 
curves and the dotted curves indicate the C x and the C y , respectively. In all cases, the side 
lengths increase in the initial stage. After such an initial stage, the side lengths reach their 
maximum values and then decrease for large ( value. When ( = 0.05, the curve does not 
show an overshoot, and the system smoothly reaches the perfect C phase. Thus, we judge 
that ( = 0.05 is the most appropriate value for our system. 



The above-mentioned dynamical SCF simulation can be performed with use of the "Sim- 
ulation Utilities for Soft and Hard Interfaces (SUSHI)" in OCTA system^. The simulation 
results reported in this article is obtained using SUSHI. 

3 Simulation Results 

We simulated the epitaxial OOT G — > C by imposing an external shear flow to an A-B 
diblock copolymer that is characterized by the parameters given in Section 2.2 using the 
technique described in the previous section. The details of the simulation procedure are 
given below. 

3.1 Initial gyroid structure and final cylindrical structure 

The initial state of the simulation is chosen as the equilibrium G structure at \N = 20. To 
generate such an equilibrium G structure, we used the following procedure. Let us denote 
the equilibrium (or steady state) side length of the unit cell of the G structure as Dg, and 
the equilibrium (steady state) spacing of the lamellar structure formed by the same block 
copolymer at xiV = 20 as Dl. The value oiD^ can easily be obtained using a one dimensional 
SCF calculation with SSO. Then, assuming an epitaxial relationship in the transitions L{10} 
— > C {10} — > G {211} at a fixed value of xN, we can obtain an approximant for Dg as follows 

D G = V6D L . (13) 

Using this value of Dq as the initial size of the simulation box, we set the SCF potential 
with the G symmetry as 

/ 2nx 2ny 2ny 2irz 2tiz 2nx\2 

V{x, y, z) = Vol cos — — sm — — + cos — — sm — + cos — sm — — , (14) 

v L>G L>G L>G Ms L>G U G ' 

where x, y, and z are the Cartesian coordinates, and the Vq is an arbitrary small coefficient 
which we assume to be 0.001 for the minor segments and -0.001 for the major segments. 
The use of the squared form on the right-hand side of eq. ()14j) originates from the fact that 
the gyroid structure in block copolymer melt is formed by double networks each with the G 
symmetry. By assigning different signs to the Vb's for major and minor segments, we can let 
the minor phase to gather inside the gyroid network while the major phase becomes rich in 
the matrix region. 

Starting from the SCF potential given by eq.(|14|). we perform a three dimensional static 
SCF calculation with SSO, which gives the equilibrium G structure. Figure 5 shows the 
optimized bicontinuous double gyroid structure obtained using the above method, where the 
parameter As in eq. (J7J) is taken as 0.2. Figure 5 shows the super cell composed of eight 
optimized conventional unit cells of the G structure. The side length of the optimized G 
unit cell is D G = 17.2. 

From this G super cell, we can extract another unit cell as shown in Figure 6(a) where X, 
Y and Z axes are chosen to be parallel to the [110], [112], and [111] directions, respectively. 
The obtained unit cell in Figure 6(a) is the minimal periodic unit cell with the Z axis oriented 
to the [111] direction of the G unit cell. The side lengths of the unit cell are \^2D G , \^6D g , 
and (\/3/2)D G , respectively, where the volume of the unit cell is three times larger than that 
of the cubic G unit cell. 



Figure 7 shows the projections of the G structure onto three different directions. The 
bicontinuously arranged rods are the domains composed of the minor A phase. Figure 7(a) 
shows the projection along the [111] direction that is the same as the left-hand side picture 
in Figure 6(a). Figures 7(b) and 7(c) show the projections along the [110] and the [112] 
directions, respectively. In Figure 7, we can see the edges of the G unit cell (tilted cube) 
drawn by dotted lines and the extracted unit cell (cuboid) drawn by solid lines. 

The self-consistent field on the X-Y plane in Figure 6(a) is used as the initial condition 
for the static two dimensional SCF calculation for the C structure at xiV = 15, where we 
assumed an epitaxial OOT G — > C. The optimized two dimensional C structure is shown in 
Figure 6(b), where we used the same scale as in Figure 6(a) for a direct comparison. 

The lengths of the vertical and horizontal axes of the two dimensional C structure shown 
in Figure 6(b) are 2.0% and 3.2% larger than those of the G structure shown in Figure 6(a), 
respectively. As the changes in the side lengths are rather minor, we expect an epitaxial 
transition for the G structure at xiV = 20 to the C structure at xiV = 15. The direction of 
the {10} plane of the cylindrical domains in Figure 6(b) coincides with that of the cylindrical 
domains in Figure 1(b). This result contradicts the standard explanation of the epitaxial 
transition G {211} — ► C {10} which was proposed in the previous experimental works and 
mean field calculations. Instead, we expect that the actual epitaxial transition should be G 
{220} ^ C {10} as shown in Figure 1(b). 

3.2 The epitaxial OOT from the G structure to the C structure 

The OOT G — > C is induced by a sudden increase in the temperature from \N = 20 to 
xN = 15, the former and the latter corresponding to the G and C phases, respectively 
This phase transition is believed to be first order and should basically be driven by the 
thermal fluctuations. An introduction of an external flow accelerates the transition 37 . We 
introduce a shear flow whose direction is oriented to the [111] direction of the G unit cell. 
The velocity field v(r) of this external shear flow is given by 

v(r) = (0,0,7(^-y)), (15) 

where y is the Cartesian coordinate along the Y axis and the C y /2 is the F-coordinate of 
the center of the system. This flow field is indicated in Figure 6(a) by the arrows. The Lees- 
Edwards boundary condition was employed in the Y direction^, and the periodic boundary 
conditions were employed in the other directions. 

For the dynamical SCF calculation, the parameters were set as follows: the criterion of 
the convergence of the segment density is A<p = 0.0005, i.e. if the difference between the two 
segment density fields at consecutive steps in the SCF iteration becomes everywhere below 
A(j), we regard the segment density field has converged. The mobility in eq. (JSJ) is set 
to Lk = 1.0, the shear rate 7 = 0.001, and At = 0.01, respectively. The parameter Q for 
the SSO is set 0.05 with which the SSO can reproduce the complete C domain in the two 
dimensional system as described in Section 2.2. With this parameter, the SSO is performed 
at every other 100 time steps. 

The temporal change of the microphase structure is shown in Figure 8. The G structure is 
deformed by the shear flow as shown in Figures 8(a) and (b). Suddenly, a grain boundary is 
generated in Figure 8(c), which is indicated by a white arrow. This grain boundary consists 
of several cylinders parallel to the [111] direction of the G unit cell and this boundary region 



separates the upper G phase and the lower G phase. The transition from the G structure to 
the C structure takes place in this lower phase as shown in Figures 8(d)-(f). The cylinders 
are tilted to the [111] direction of the G unit cell as shown in the side view of Figure 8(f). 

The tilting of the cylinders is caused by the constant shear flow because a steady shear 
flow is composed of two contributions, a uniaxial extension and a rotation This rota- 
tional contribution tilts the cylinders. Such a tilting is suppressed in experiments by using 
oscillatory shear flows. 

Three fold junctions in the upper G phase shown in Figure 8 (f) are stable and are 
migrated by the shear flow. Even in the late stage (t— 10000), we cannot obtain the final 
equilibrium C structure. This result suggests that boundary between the upper G phase and 
the lower C phase is stable and the separated G and C phases coexist stably. Actually, a 
clear boundary between the G and C grains is observed by a polarized optical microscopy 
in a polymer solution an d the long-lived coexistence between the G and the C phases is 
experimentally observed in a block copolymer^. 

The changes in the side lengths of the simulation box are shown in Figure 9. The side 
lengths in the X and Z directions are almost constant, which means that the epitaxial 
condition is satisfied. On the other hand, the side length in the Y direction increases with 
time. The reason of the increase is explained below. 

To check the effect of the SSO method, we carried out the same dynamical SCF simula- 
tion but without the SSO method. The time evolution of the domain morphology is shown 
in Figure 10. In this shown in Figures 10(a)-(d), the OOT occurs at the center of the 

system where the G structure transforms into the C structure. The stable grain boundary 
as shown in Figure 8 is, however, not observed in this case without SSO. The cylindrical 
domains are also tilted to the [111] direction of the G unit cell as is shown in the side view 
of Figure 10(c). After such a transient state, the system reaches the complete C structure. 
A characteristic phenomenon is observed near the center of the system where the cylindrical 
domains reconnect as shown in Figures 8(d)-(f). Such reconnections continue steadily for 
certain time duration. This reconnection phenomenon means that the system is in a dynam- 
ical steady state where the energy injected by the shear flow into the system is released by 
the energy dissipation accompanied by the periodic reconnections of the cylindrical domains. 

In order to check the stability of the systems, we show in Figure 11 the time evolution of 
the free energy density during the phase transition. The free energy density in the case with 
the SSO shows a moderate change compared to that in the case without the SSO, which 
shows an oscillation synchronized to the reconnections of the cylindrical domains. Such a 
periodic change in the free energy density is observed after t = 3000 when the system reaches 
the almost perfect C structure without defects. In order to obtain such a perfect C phase, the 
system goes over energy barriers by the driving force of the shear flow. In the case without 
SSO, however, the energy barrier should be much higher than the case with SSO because 
the condition of the constant system size imposes a sever restrictions on the reconnection of 
the cylindrical domains. In the case with SSO, such a restriction is avoided by the increase 
in the side length of the simulation box in the Y direction as shown in Figure 9. 

We also tried simulations under a shear flow whose velocity gradient is set parallel to the 
X direction. In this case, the free energy of the system increases slightly but the nucleation 
and growth of cylindrical domains can not be observed even in the late stage £=6000 either 
with SSO or without SSO. The epitaxial condition for the OOT is expected to be more pre- 
cise for this direction of the velocity gradient than the case with the velocity gradient in the 
Y-direction. This is because the periodicity of the G structure in X direction matches the 



periodicity of the C structure than that in the Y direction, the former promoting the gener- 
ation of the C structure. Contrary to this expectation, three fold junctions perpendicularly 
oriented to the [111] direction of the G unit cell continue to disconnect and reconnect due to 
the shear flow. Figure 12 shows this phenomenon. The circles in the Figure 12 indicate the 
three fold junctions where the disconnections and the reconnections take place. Figure 12(a) 
shows the structure after the disconnections, where we can observe the remains of three fold 
junctions indicated by the circles. Figure 12(b) shows the structure after the reconnections, 
where the three fold junctions regenerated. This result indicates that the G structure has 
different stabilities to different directions of the shear velocity gradient. 

The reason of this different stability is explained by using Figure 13. Figure 13 shows 
three fold junctions perpendicularly oriented to the [111] direction of the G unit cell with 
different rotational angles to the direction of the shear gradient. Figure 13(a) shows a three 
fold junction under a shear flow with its velocity gradient in the Y direction. The three 
domains extending from the center of the three fold junction are subjected to different shear 
flow velocities, i.e. the three domains do not move with the same velocity v y . Thus, the 
three domains are elongated to the different directions with different rates, and the elongation 
finally makes the three fold junction disconnected. On the other hand, in the case of a three 
fold junction under a shear flow with its velocity gradient in the X direction as shown 
in Figure 13(b), two domains have the same velocity v x . In this situation, the three fold 
junction is elongated to the positive direction of the X-axis. Even after the disconnection, 
the separated domains keep closer and will be reconnected easily to form the three fold 
junction structure as shown in Figure 12. 

4 Discussions 

Our simulation showed that the epitaxial OOT G — > C takes place in the [111] direction 
of the G unit cell and that the epitaxial relation for the G {220} — > C {10} transition is 
achieved. The transition does not occur uniformly as shown in Figures 8(d)-(f) and Figure 
10(c), where the C domain nucleates and grows. 

Most of the experiments have reported the epitaxial OOT C {10} <-> G {211} which 
disagrees with our simulation result. A possible reason of this discrepancy is as follows. The 
first diffraction peak from a G structure is the peak from {211} and the intensity of this 
peak is stronger than that of {220} peak (secondary peak). Moreover, the positions of the 
peaks from G {220}, G {211} and C {10} are so close that it is not easy to judge which of 
the peaks from G {220} and G {211} epitaxially matches with the peak from C {10}. We 
calculated a three dimensional scattering function of the optimized G structure obtained in 
the simulation, and confirmed that the {211} spots are dominant and their intensities are 
about four times larger than those of the {220} spots. 

If the G {211} — > C {10} is realized under a shear flow with the velocity gradient in the 
Y direction, the planes composed of cylinders are directed in parallel to the sheared plane, 
i.e. the XZ plane. Thus, the friction generated by the reconnecting domains to the shear 
flow is expected to be smaller than that for the G {220} — > C {10} case where the cylinder 
planes are perpendicular to the sheared plane. In our simulations, however, the system 
prefers the pathway as G {220} — > C {10}. Therefore, we conclude that the direction of 
the velocity gradient is not an important factor in determining the direction of the C planes 
for the epitaxial transition. On the other hand, we confirmed that matching between the 



lattice constants is more important. As is shown in Figure 6(b), the origin of the selection 
of the generated C structure from the G [111] plane (Figure 6(a)) is the matching between 
the lattice constants. That is, the system prefers the kinetic pathway that minimizes the 
free energy of the system by matching the lattice constants, which leads to the G {220} — > 
C {10} transition. 

Previous theoretical studies have also supported the OOT G {211} — > C {10}. These 
studies relied on the reciprocal space representations. However most of the experiments have 
been done under a condition with a shear flow and a temperature change. We succeeded in 
reproducing such experimental conditions in our simulation. Using this simulation, we could 
reproduce the correct kinetic pathway of the epitaxial OOT, i.e. the nucleation and growth 
process of the C domains. 

We found the difference in the stability of the G domains to the shear gradient direction 
due to the different velocities of the shear flow imposed on the three domains meeting at a 
three fold junction as shown in Figures 12 and 13. The G structure is stable under the shear 
flow with the velocity gradient in the X direction. There has been no answer to the question 
why the complex G phase with three dimensional bicontinuous structure is generated under 
a shear flow. Our result demonstrates that the G structure is actually stable under a shear 
flow. 

Although the detail of the transition process is complex, we understand that the three 
fold junctions with domains perpendicular to the [111] direction of the G unit cell do not 
play an important role in the transformation from G to C. Three fold junctions are simply 
disconnected and vanish during the phase transition. This observation does not agree with 
the model of the epitaxial transition proposed by Matsen, where a three fold junction is 
connected to one of the nearest neighbor three fold junctions to form a five fold junction. 
In our observation, three fold junctions are stable and they are not connected to any other 
junctions. 

Here, we propose a model of the kinetic pathway. Figure 14(a) shows a projection of the G 
unit cell along the [111] direction, where the bold triangles are the projections of consecutive 
three domains. Figure 14(b) is the same structure as Figure 14(a) observed from a different 
direction. We can confirm that the triangles in Figure 14(a) are formed by consecutive three 
domains (shown in black) connected by three fold junctions. When a shear flow is imposed, 
these black domains in Figure 14(b) are elongated and form cylinders as shown in Figure 
14(c). These cylinders are rearranged to form a hexagonally packed cylindrical structure 
whose lattice spacings satisfy the epitaxial relations G {220} — > C {10} as shown in Figure 
14(d). This model of the kinetic pathway can be verified in Figures 8(d)-(e) and in Figures 
10(b)-(c). 

5 Conclusion 

The epitaxial OOT G — > C was studied using the real space dynamical SCF technique with 
the SSO method. With such an SSO method, we succeeded in reproducing the realistic 
kinetic pathway of the first order phase transition of G — > C. On the other hand, in the 
absence of the SSO, we found that the kinetic pathway is very different from what we 
observe with SSO. We also found that the G structure shows different responses to different 
directions of the velocity gradient of the shear flow. 

Using this technique, we studied the kinetic pathway of the G — > C transition induced 



by a shear flow in the [111] direction of the unit cell of the G structure. We observed the 
following kinetic pathway: The G domains perpendicularly oriented to the [111] direction 
of the G unit cell do not contribute to the formation of the cylindrical domains. They are 
disconnected and vanish during the transition. On the other hand, the other G domains are 
elongated by the shear flow and transform into the cylindrical domains. Such deformations 
occur locally, and the cylindrical domains are rearranged to form a hexagonally close-packed 
C structure. 

The most important result of our simulations with SSO is that we can observe a nucleation 
and growth of the C phase in the matrix of the G phase, i.e. the nature of the first order 
phase transition, which was not observed in the previous simulations done in the Fourier 
space. Under a steady shear flow, we observed that the G domains around the nucleus of the 
C phase deform and gradually join the C phase. We also observed that the domain spacing 
satisfies the epitaxial relationship G {220} — > C {10} as was proposed in the experimental 
work — . We could not observe the transformation process of the domains from three fold 
junctions to five fold junctions as was previously proposed^. 

We found that the dynamical SCF theory with the SSO method in real space is very 
useful and reliable to trace the OOTs and ODTs between the microdomain structures of 
block copolymer melts. 
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Figure captions 



Figure 1. A projection of the G unit cell structure observed from the [111] direction ^31. 
The circles indicate the positions of the cylindrical domains in the epitaxial transition (a) G 
^{211} —* C G?{ 10 } and (b) G d{ 2 20} - > C <i[ 10 }, respectively. 

Figure 2. A comparison of the domain morphologies of an A-B diblock copolymer ob- 
tained with the two dimensional dynamical SCF simulations with SSO. The simulations 
were started form a D phase without the external flow. The model parameters are shown in 
the text. The values of the SSO parameter ( are (a) ( = 0.001 and (b) ( = 0.05, respectively. 

Figure 3. The time evolutions of the free energy per chain for the C structure for sev- 
eral ( values are shown. The dotted curve is the result of the reference simulation with 

C = o. 

Figure 4. The time evolutions of the side lengths of the simulation box for several ( values 
are shown. The solid and dotted curves show C x and C y , respectively. Both side lengths are 
normalized using their initial values C{ = 32 as £i/£i . 

Figure 5. The isosurfaces of the super cells (2 x 2 x 2) of a bicontinous double gyroid struc- 
ture for xN — 20, (f) = 0.9 are shown. The single unit cell is calculated using 32 x 32 x 32 
meshes. 

Figure 6. A comparison of the G structure and the two dimensional C structure obtained by 
the static SCF simulations with SSO is shown. These two structures indicate an epitaxial 
OOT G — > C in the [111] direction, (a) A view of the isosurfaces of the minimal periodic 
structure of the G unit cell along the [111] direction is shown, The parameter are chosen as 
xN = 20, and = 0.75 (left-hand side). The velocity distribution of the added shear flow 
in the [111] direction (right-hand side), (b) The two dimensional C structure obtained as an 
equilibrium state starting from the domain structure on the cross section on a {111} plane 
of Figure 6 (a) with xN = 15. 

Figure 7. Projections of the real space G structure onto different directions are shown, 
where the bicontinuously arranged rods correspond to the minor A phase of the G structure. 
Shown are the projections along (a) the [111] direction, (b) the [110] direction, and (c) the 
[1 12] direction, respectively. In (a), two different unit cells are shown: one is the conventional 
cubic unit cell drawn by dotted lines and the other is the parallelepiped unit cell which has a 
periodicity in the [111] direction drawn in solid lines. The relation between the side lengths 
of these two unit cells are discussed in the text. 

Figure 8. The time evolution of the domain in the epitaxial OOT G — > C simulated with 
SSO is shown. A shear flow is imposed in the [111] direction of the G unit cell with 7 = 0.01. 
The directions of the shear flow are shown by black arrows. The view point of the graphics 
is set so that one can verify the growth of the C structure. The snapshot figures are taken 
at times t = (a) 100, (b) 960, (c) 1960, (d) 3960, (e) 4960, and (f) 5960, respectively. For 
the cases (c) and (f), we also show the side views, where the white arrow indicates the stable 
grain boundary where cylindrical domains parallel to the shear direction ([111] direction) is 



generated. 



Figure 9. The time evolutions of the side lengths of the simulation box during the phase 
transition are shown. 

Figure 10. Similar to Figure 8 but without SSO. 

Figure 11. The time evolutions of the free energy per polymer chain with SSO and without 
SSO during the phase transition are shown. In both of these cases, the initial value of the 
free energy is set to be zero. 

Figure 12. The time evolution of the G domain under a shear flow with its velocity gra- 
dient in the X direction and with 7 = 0.01. The circles indicate the three fold junction 
perpendicularly oriented to the [111] direction of the G structure, (a) Three fold junctions 
are disconnected at £=3600. (b) Three fold junctions are regenerated at t=4200. 

Figure 13. Three fold junctions are migrated to different directions by the shear flow, (a) 
Under a shear flow with velocity gradient in the ^-direction and (b) a shear flow with veloc- 
ity gradient in the X-direction. In these two cases, the domains migrate in different manners. 

Figure 14. The model of the kinetic pathway of the epitaxial OOT G — > C is illustrated, 
(a) Shown is the view of the G unit cell along the [111] direction, where the black domains 
transform to the C domains, (b) Another view of the same G unit cell as in (a) from a 
different direction is shown, (c) An image of the transformation from the elongated domains 
to the cylinders is shown, (d) Positions of the rearranged cylinders with hexagonal order are 
shown, where the lattice spacings satisfy the epitaxial relations G <i{22o} C d{ W y. 
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